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The diffusive properties in velocity fields whose small scales are parameterized by non 5-correlated 
noise is investigated using multiscale technique. The analytical expression of the eddy diffusivity 
tensor is found for a 2D steady shear flow and it is an increasing function of the characteristic noise 
decorrelation time t. In order to study a generic flow v, a small-r expansion is performed and 
the flrst correction 0{t) to the effective diffusion coefficients is evaluated. This is done using two 
' different approaches and it results that at the order r the problem with a colored noise is equivalent 

^\ , to the white-in-time case provided by a renormalization of the velocity field v ^ v depending on 

0^ ' T- Two examples of 2D closed-streamlines velocity field are considered and in both the cases an 

' enhancement of the diffusion is found. 
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CN ; I. INTRODUCTION 

^ i 

^ , The problem of diffusion in a given velocity field has both theoretical and practical relevance in many different 
[ — ' fields of science and engineering as, e.g., transport processes in chemical engineering and combustion studies [Q. 
^\ [ The tracers transport, in particular the evolution of their concentration, plays an important role in many aspects of 
■ geophysics. For the oceanic flows, satellite data indicate that the mesoscale features, like eddies and cold filaments, 
advect temperature and nutrients over spatial and temporal scales longer than those of the geostrophic turbulence. 
The diffusion enhancement by a given velocity field has attracted a lot of works in the last years. In particular the 
role of the velocity field properties has been largely investigated while the effects of small scales parameterization are 
not understood. 

' In this paper we will focus on the effects of a finite noise correlation time. This problem is relevant in studying the 
" transport in the ocean since in this system the noise term comes from unresolved velocity scales which are correlated 
I in time. 

' O In Section 2, by using the multiscale technique, we study the diffusion properties of the model proposed in Ref. 

1^ for transport in the upper mesoscale ocean. The transport is described by a Langevin equation with a Gaussian 
^ ' colored noise in time. 

. . , The aim is to understand whether a finite noise correlation time r enhances or depresses the dispersion process in 
^ ' a given velocity field v{x, t) with respect to the delta-correlated case (t = 0). 

k>( [ Exploiting the scale separation in the dynamics we derive, using the multiscale technique |^, an effective diffusive 
j_] ' equation for the macrodynamics, the calculation of the effective diffusivity second-order tensor is reduced to the 
d I solution of one auxiliary partial differential equation |^ , Q , |8| . 

In Section 3 we consider a shear ffow, in this case the diffusion coefficient increases with r. The solution of the 
auxiliary equation is, in general, quite difficult, therefore, to investigate the role of the finite t in Section 4 we perform 
a small-r expansion. An alternative method is presented in the Appendix A. 

In Section 5 we study the case of two closed-streamlines fields that mimics the transport in the Rayleigh-Benard 
system: the quasi-two-dimensional flow studied by Shraiman in [|l^ and the AB flow. In both the cases the presence 
of a small correlation time enhances the diffusion process. 
Conclusions are reserved for the flnal Section 6. 

II. EFFECTIVE DIFFUSION EQUATION FOR GAUSSIAN COLORED NOISE 

We consider large scale and long time dynamics of the model proposed in ||^ and already studied in ||^ for the 
transport of a fluid particle in the upper mesoscale ocean: 

^x^v{x,t) + s{t) (1) 



o 

0^ 



1 



where v is a cZ-dimensional incompressible velocity field (V ■ v = 0), for simplicity, periodic both in space and in time 
and s is a Gaussian random variable of zero mean and correlation function 



{s,(t)s,{t'))^-6,,e-'^. (2) 

T 

The term v{x,t) represents the part of the velocity field that one is able to resolve, i.e., the larger scale mean flow, 
whereas s{t) represents the part of the velocity field containing the subgridscale flow processes, e.g. the small-scale 
turbulence. The plausibility of such a description is discussed in [0, ||Tl[] . 

In the hmit r — > 0, resulting e^l*~*'l/'^/T S{t — t'), the (§) reproduces the widely studied delta-correlated case 

(s,it) sj{t')) =a^ 6^,Sit~t') (3) 

the diffusive properties of which we would like to compare with the r-correlated noise case. 

To study the dispersion of tracers evolving according to eqs. (|I|) and (||) on large scales and long times we use the 
multiscale technique. This is a powerful mathematical method, also known as homogenization, for studying transport 
processes on time and spatial scales much larger than those of the velocity field v. It has been already applied to the 
delta-correlated case [|j and it has been shown that the motion on large time and spatial scales is diffusive and it is 
described by an effective diffusion tensor which takes into account the effects of the advecting velocity on the bare 
diffusion coefficient a^. 

To apply this method to the case of Gaussian colored noise, we first write eqs. (0) and (^) into a Markovian process 
by enlarging the state space considering s{t) as a variable evolving according to the Langevin equation: 

d 1 , ^ , 

-s = —s + w{t) 4 
at T 



where now the noise w{t) is a white noise with correlation functions 

, 2 



{w,it) w,{t')) = 2[^y 5.,,5{t-t'). (5) 



We have now a two- variable {x, s) Markovian process whose associated Fokkcr-Planck equation can be easily obtained. 
Indeed introducing 



the equations (0) and become 



The associated Fokkcr-Planck equation is 



^Y = V{Y,t)+A-W{t). (7) 



= (7) ' dss^Q -{v + s)-d^Q + ^ds- [s e) (8) 

where <d — Q{x^s,t) denotes the probability density. 

The doubling of the space dimension is the price to pay for having a Fokker-Planck equation. In the Appendix A 
we discuss a different approach to the problem which does not double the dimension of the space, but leads in general 
to a non-Markovian master equation. 

We can now apply the multiscale technique. Following in addition to the fast variables x and t we introduce the 
slow variables defined as X = ex and T = e^t where e <C 1 is the parameter controlling the separation between the 
small scales related to the velocity field v and the large scale related to the Q variation. The two sets of variables are 
considered independent and so we have to make the substitution 

d^^d^ + tdx ; dt^dt + e^dr- (9) 

The solution of the Fokker-Planck equation (||) is sought as a perturbative series 

e(a;, t, X, T, s) = G^^^ + eO^^) + e^e^^) + ... (10) 
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where the functions 8*^"-' depend on both fast and slow variables. By inserting ^ and (|o|) into the Fokker-Planck 
equation (|^), equating terms of equal powers in e and choosing the solutions which have the same periodicities as the 
velocity field, we obtain a hierarchy of equations the first three of which are: 

X'e(")=0 (11) 
Pe^i) = -(t; + s)-9xe(°) (12) 

i?e(2) = -(t; + s).9xe(i)-9Te(°) (13) 

where the operator V is defined as 

V = dt + {v + s)-d^-^dsis )-(^)'as/. (14) 

In order to solve eq. (^ij) we make use of scale separation and we write the solution 0^^^ as the sum of two terms: 
(0'^*'^)(X, r, s) depending on the slow variables and Q^'^\x,t, s) depending on the fast variables. Here and in the 
following the ( • ) indicates the average over the fast variables. 
Equation ( pi] ) then splits into the two equations 

V e(°) = (15) 

^Ydj{ei")) + ^d, (s(e("))) = o. (16) 



One can show that the solution 8*^°-' will relax to a constant with respect to fast variables, so we can simply take 

e(o)(a;,i,X,T,s) = (e(°))(X,T,s) (17) 

and write the solution of eq. ( |l6|) as 

(e("))-Q(X,T)P(,s) (18) 

where 'P(s) is defined as 

= , (19) 

(zttct^/t) 2 

with d the dimension of the a;-space. 

By using ( p^ we see that eq. ( |l2| ) can be written as 

=/(a;,i,s).G(X,T) (20) 

with 

f{x,t,s)=-{v + s)V{s) , G{X,T) = dxQ{X,T) (21) 

with solution 

e(i)(a;, t, X, T, s) = x{x, t, s) ■ G{X, T). (22) 
The vector field x is called the auxiliary field and it solves the auxiliary equation 

'Dx = f- (23) 
Finally by averaging eq. ( [l^ over the fast variables , ( ■ ) , and integrating over s , ( • ) , we obtain 



dx{{v + s) 9(1)) =-aT(e(o)) = -dTQ{X,T) (24) 

which, using eq. (|2^), becomes 
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OtQ = -{{v, + s,) x^) ai^x, Q = Df^ dhx, Q- (25) 

This is the diffusion equation describing the large-scale dynamics, i.e., the dynamics in the slow variables. The 
effective eddy diffusivity tensor Df^ is given by 

Df,=-\ (((^^TT^^yxTT + {{vj + X. )) ■ (26) 

From the auxiliary equation (^3|) one can show that the Df^ is positive definite. Indeed, if we consider the «-th and 
the j-th component of (^3|), multiply by Xj ^-nd Xi respectively, sum the two terms and average the result over the 
periodicities, ( • ), and integrate over the random variable s we end up with 

= (7)'/ d'^s i^s{x^)) (dAXj)) V{s) > 0. (27) 
This result can be extended to non-periodic velocity field following the prescriptions in m . 



III. A SOLVABLE CASE: THE STATIONARY SHEAR FLOW 



The resolution of the auxiliary field equation for a generic v is not an easy task. Therefore not trivial solvable 
cases are useful to understand the properties of the solution. In particular the auxiliary equation can be resolved for 
parallel flows, which in two dimensions have the form 



v{x,y) = iv{y),0) 



(28) 



where v{y) is an arbitrary function of y. Note that these flows automatically satisfy V • 7; = 0. To evaluate the 
effective diffusion coefficients we first write the solution of th auxiliary equation as 



where Q is the Green function of the operator T): 

Vg = S{x- x') 5{t - t') 5{s - s'). 

Inserting (E9[) into eq. (Eq) we have: 



Now we note that the Green function Q can be written as 



{5{x - x{t; x\ <')) '5(s - s{t; s\ t')))i 



(29) 



(30) 



(31) 



(32) 



where the average is over the realizations of the white noise ii;, x[t\ a;', t') and s(t; s', t') are the solutions of eqs. (|^), 
(^ and (^Shwith initial condition x' = x{t' \ a;', i') and s' = sit' \ s' , t'). For the velocity field ( p8| ) the solutions of the 
equations^ and (^ can be written as 

x{t) ^x' + dti v{y{ti)) + //, dh si(ii) 

y{t) = y' + J', dh S2{ti) (33) 
s{t) = s'e-(*-*')/^ + dti ^(;(^l)e-(*-*l)/^ 

Inserting and ([3^) into (p^), after some straightforward straightforward algebra one obtains: 



27r 



Df^=a^ + — dk \ d{k) p lim / dt' exp (-a^ (t - t') - r (l - e 



(34) 



and 



1^22 = O- 



(35) 
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where v{k) is the Fourier transform of v{y). The same result can be obtained directly from the definition 

Dfj = lim hix^it) ~ (x,)) (xjit) - (xj))) 

^ t — >oo zt 

using (^3|). Now, because of the inequality 

exp 



> 1 



(36) 



(37) 



one has 



i?fi(r) > ^2 + ^|dfc I vik) f \\mj^ 

V{k) 12 



dt' e-^" 



27r 



dk- 



(t2 fc2 



(38) 



Therefore for a stationary parallel flow a colored noise produces an enhancement of the dispersion. Similar equations 
can be obtained for a time dependent shear flow. However in this case it is not simple to see the sign of the correction. 
The results will be reported elsewhere. It is trivial to show that the results in this section hold also for a 3-d shear 
flow : 



v(a;,j/,z) ^ (i;(y,z),0,0) 



(39) 



IV. EDDY DIFFUSIVITY FOR SHORT NOISE CORRELATION TIME 



By using multiscale technique, the calculation of the eddy diffusivities has been reduced to the solution of the 
auxiliary equation (23). Numerical methods are generally needed to solve it but to do so we have to work in a 
2d-dimensional space. In general, this is not feasible so, to get more insight the generic v case we study the small 
r case and expand the auxiliary field x in a power series of r. A typical time of the physical system which can be 
compared to r is = A/ {v^)^/"^ i.e. the average time it takes a particle to travel a characteristic length A. 

Taking 



'T Si 
T 



fc=0 



we obtain the following expression for the eddy diffusivity tensor 



= 

2 27ra2 



. k 



and for the auxiliary equation: 



E 



[rt [dt+v- a,) + a • 0, - (a„(a ) - a^O, 



2^2 ^ 

aa J 



.(fe) 



(40) 



(41) 



From the expression ( }40D we see that in order to determine the correction of order r to Df< , we need the quantities 
{X^^°^), Xl^°\ T^^'\ {X^^^^), Xl^^\ (x.^'^)- The fields x.('=) obey to the equations 



(0) 



(0) 



a • dxX 



(1) 



+ Vie 2a 



/i > 3 



(42) 

(43) 

(44) 
(45) 
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where the operator Oa is defined as 

The solutions of cqs. (^), ( ^ ) and (|4^ ) can be written in the form 

2 



(2) 



{biia\ + b2ia\ + 6310102 + 64101 + 65ia2 + bm 



the coefficients a^j and 6.^ are functions of x and i while for the field Xi"^^ it results that 



(46) 

(47) 
(48) 
(49) 

(50) 



By inserting expressions (f47|), ( f4S| ) and (^) into eqs. (|42|), (|43|), ( |44[ ) and (|45|) we can determine the coefficients a^- 
for I = 1, 2 and 5^ for i = 1, ••, 5 and then by integrating over a the equations for Xj^'^\ Xj^^^ s-^id Xj''^'' respectively 
we finally have the equations for the remaining functions xf^ i ct3i ^ind b^i 



0.,t x/°^ = 

Oxt 031 = 

3 2 



^\^+^\^JX^ 



2 cr-9ji;„ a^^Xi 



where 

Oxt = dt+v-d- a'd 
The D,^ coefficients at the first order in r read 



2n2 



(51) 
(52) 

(53) 



(54) 



- ^ Xf > + {V, xf ))) 



(55) 



We note that instead of the 2(i-dimensional equation ( |26| ) we have now a system of two d-dimensional equations ( [5l| ) 
and (|5^ ) without the random variable s. Of course this is numerically much more convenient. 
We note that by defining new velocity and auxiliary fields as 



cr T ^2 

V = V O V 

2 



and negletting terms of ©(t^) equations (pT]), ( [S^ ) and ( |55| ) can be written as 

(i9t + i>9 - (7^5^) Xi = ^'i 

and 



(56) 
(57) 

(58) 
(59) 



formally equivalent to the Gaussian white noise result. 

In the Appendix A we use a different method to obtain the expression of the eddy diffusivity tensor for small r 
up to 0{t^). By starting from the Master equation associated with the Langevin equation (0) and using a small t 
expansion we derive the Fokker-Planck equation 



dte{x,t) = -d, [v,{x,t) e{x,t)] + dl, [v,,{x,t) Q{x,t)] 



(60) 
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with 



Now, by applying multiscale technique to (|60[) we end up with the following equations 



w 



(0)\ 



(0), 



„(1)\ 



„(1) 



(61) 

(62) 
(63) 

(64) 



which differ from the previous ones. 

This result poses the question about the validity of the two sets of equations and hence of the expansions. In order 
to answer to this question we firstly note that the two sets can be considered as a particular case of the generalized 
equations 



O^twi'^ = a' ((p - q) + {p-q) d.wf^ +pduv, d^wf'^) 



T 

~ 2 



obtained starting from the Master equation 



where 



dte{x,t) = -v^{x,t) ds{x,t) + df^ [v,j{x,t)e{x,t)] 



Vi = Vi + q (T^rd^ Vi 
Vij = o-^ (Sij + ^ {di Vj + dj Vi) 



(65) 
(66) 

(67) 

(68) 

(69) 
(70) 

and by applying the multiscale technique. 

The role of the two free parameters p and q is of generalizing the Master equation (|6^) being our feeling that 
there exists a family of different microscopic process specified by v and Vij that correspond to the same macroscopic 
diffusive process specified by Dfj. 

The first set consisting of the equations 
second one consisting of the equations f|62| 



( p3| ) and ( |5q ) corresponds to g — 1/2 and p — 2, whereas the 

and (^4|) corresponds to q = and p — I. Now a closer analysis of 

the diffusion coefficients ( |70| ) reveals that it may take negative values, introducing unphysical singularities into the 
problem. Moreover the multiscale technique requires a Df^ definite positive. This means that ( |68| ) do not constitute 
genuine Markovian process with a well defined corresponding Langevin equation driven by a Gaussian white noise. 
In general this expansion do not converge uniformly in x and the range of validity is restricted to r <C 1, r/a^ ^ 1 
(in dimensionless unit), thus they are asymptotic estimate for t —f 0. In this range of equations validity we do not 
have any a priori arguments for choosing one of the two sets of equations and it is for this reason that we expect to 
obtain the same value of the eddy diffusivity tensor using the two sets. This is indeed the case for the parallel shear 
flows (see Appendix B). This result can be thought as a first indication of the equivalence (with respect to the value 
of the diffusivity) of the generalized equations. 

In general we do not expect that all possible choices of q and p give the same Dfj but we can think to select the 
class of equivalent process by applying the generalized equations to the parallel shear flow and imposing Dfj equal to 
the known expression ( ^8|) independent on both q and p. This calculation is reported in Appendix B and it ends up 
in the following condition 



p^2q+l 



(71) 
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All these considerations suggest that for short t there exists a class of equivalent equations depending on the 
parameter p that lead to the same eddy diffusivity tensor. 

If this is the case, among the all possible choices of q andp consistent with (O), we can choose: p — Q and q = —1/2. 



In this case 'Dij — a Sij and (60) reduces to a truly Markovian process with the associated Langevin equation 



j^x^vix,t)+m (72) 

where ^ is a Gaussian white noise. This is the microscopic Markovian process which approximates the long time 
and large space transport properties of a colored noise process with short noise correlation time. In other words to 
study the diffusion properties of (|l|)-(||) for small r we can replace the original colored noise process with the process 
described by ([72[). This will give the correct diffusion coefficients up to C(t^). 

We have checked numerically that for the AB flow the eddy diffusivity tensor assumes the same value for the three 
different choices of the parameters q and p consistent with ([7l|): p ^ 2 and q = 1/2, p = I and q ^ 0, p = and 
q = —1/2. This make us confident with our conclusions. 

V. FLOWS WITH CLOSED STREAMLINES 

We apply now our analysis to two models for the Rayleigh-Benard steady convection: the first one consists of 
an horizontal extent of convection cells much larger than its height so that the flow can be considered quasi-two- 
dimensional; the second one is the two-dimensional AB flow made of a structure periodically repeated in the space. 

A. A quasi-two-dimensional flow 

We consider the flow discussed by Shraiman . This is described by the stream function 

, , vL . f nk \ , / TT \ 
V[x,y) = —rsm[ —X ] sm[-y] (73) 



nk \ L J \L 

with V being the characteristic velocity, L the height of the cell (y e [0, L]) and L/k the x-periodicity of the roll 
pattern. The top and the bottom plates of the cell are assumed impermeable for the passive scalar so that the 
appropriate boundary conditions for the tracers density function G are dyQ \y=a.L= 0. The streamlines of the flow 
are illustrated in Fig. (|^). Using Fokker-Planck equation 

dte^-v-de + a^^d^e, (74) 

for large Peclet number |^ Pe = vL/a'^, k = 2 and r = the eddy diffusivity coefficient £'f]^(0) has been calculated in 
ItI and it is 



VTT V TT 

According to the results of the previous section, the colored noise case with small t is described up to order O(t^) 
by the same Fokker-Planck equation provided the velocity field is renormalized as eq. ( |69| ) with q — —1/2. Taking 
into account that Vx = —dy^, Vy = dxip, and using ([73|), we have 

5V = -(7r/i)^ (l + fc2) V- 
Therefore up to order 0{t^) the diffusion is described by 



'^The case studied by Shraiman corresponds to the delta-correlated case { r — 
where is now the molecular diffusivity of the considered system. 
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^{x, y) = (l + ^ {j^y (l + ^(^^ y) = c (a^, r, fc) y). 



(76) 



Since c does not depend on x and i, we can repeat the calculation of Shraiman and obtain under the same conditions 
of (^) the expression for the eddy diffusivity coefficient 



^2 



5 cr^r / TT 



2?fi(r) = -^^/^=i5fi(0)^(^l 

= Dfm + (77) 
We then conclude that, in this small T enhances the diffusion coefficient. 



The same result can be deduced from the multiscale equations (58) and (^9|): in fact because of the structure of 
these two equations it is not difficult to show that if we change only the module of the velocity field {v = cv) and we 

Qow the explicit form of Df^ = f{v) as a. function of v we have Df^ 
function f{v) is given by eq. (^5|) from which (^^ follows. 



know the explicit form of Df< = f{v) as a function of v we have Df- — f{v). For large Peclet number and t = the 



B. The AB flow 



The AB flow is given by the velocity field 



v{x,y) = (i?cos(y), v4cos(a;)). 



(78) 



For A = B = 1 the streamlines form a closed periodically repeated structure made of four cells as shown in Fig. (|^). 
We expect that the diffusive behavior of such a system is similar to the previous case, in fact we know that for 
small Peclet number Pe the eddy diffusivity tensor is proportional to V Pe like in the quasi-two-dimensional case. 
In the Figure @j the behavior of the quantity A = [D'^{t) - D^Q)]/[D^) r] versus is shown. The correction 



A has been calculated by integrating numerically the equations ( pi]) and ( |53| ) and evaluating the quantity 



DS(0) 



(0)\ 



{vi b, 



'61/ 



A 



(79) 



for different values of a^. The equations are solved by using a pseudo-spectral method ||T^ in the basic periodicity 
cell with a grid mesh of 64 x 64 points. De-aliasing has been obtained by a proper circular truncation which ensures 
better isotropy of numerical treatment. 

It is evident that also in this case the numerical results follow a linear behavior with a positive angular coefficient 
and so we can conclude that the introduction of the colored noise leads to an enhancement of the diffusion. In 
particular we can see that the numerical results follow very well the line A = cr^/4. This is not surprising because we 
know that for large Pe 



therefore using the same arguments of the previous section we can deduce that 



(80) 



(81) 



in a very good agreement with the numerical results. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have studied the transport properties in velocity fields whose small scales are parameterized by 
Gaussian colored noise. We analyzed in particular the effects of a finite noise correlation time r on the diffusive 
properties for large time and spatial scales. In this limit, using the multiscale technique, we derive the diffusion 
equation ( |25| ) and the associated effective diffusion tensor. The latter is obtained, once the velocity field v is given, by 
the solution of an auxiliary equation, see (|^ and (p6|), of the same structure of the original Fokker-Planck equation 
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| [19[ . The former is, however, an exact result for the diffusive regime vahd for very long times, thus avoiding all finite 
time effects of the Fokker-Planck equation or the associated Langevin equation 

The auxiliary equation cannot be solved for a generic velocity field, nevertheless there are nontrivial flows for which 
the solution can be found. This is the case for the steady parallel flow for which the effective diffusion coefficient is 
an increasing function of r. To study in more details the effects of a small noise correlation time for a generic v we 
have performed a small-r expansion and evaluated the first correction 0(t) to the effective diffusion coefficient. This 
is done by using two different approach. We find that to order 0{t) there exist a one-parameter family of flows with 
the same diffusion properties. This invariance can be used to pick up the most convenient microscopic dynamics, 
from both analytical and numerical porpoises. We apply the small-r results to two two-dimensional model flows with 
closed streamlines. In both the cases we find an enhancement of the diffusion. 

The enhancement of the diffusion for small r has been interpreted in p8| in terms of interference mechanism 
between turbulent and molecular diffusion. The colored noise makes the diffusion particles forgotten of their previous 
positions less rapidly than in the white noise case, thus the Lagrangian correlation time increases and so does the 
eddy-diffusivity. 

The study of the problem for t not small is object of current work. 
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APPENDIX A: THE MASTER EQUATION FOR COLORED NOISE AND SMALL r 

In this appendix we derive the master equation for the probability density Q{x,t) in the limit of small r for the 
process described by the Langevin equation (|^) and If the random variable s is a Gaussian white noise of zero 
mean, 8(a;, t) satisfies the Fokker-Planck equation. We address here the case of a colored noise. Now the process x is 
non-Markovian and no exact simple equation for Q is known. Let us consider an Ornstein-Uhlenbeck process s, that 
is a zero mean Gaussian process with correlations [c.f.r eq. (0)] 



(Al) 



where r is the correlation time. The probability density is given by 

e(a;,i) = {Six{t)-x)) 



(A2) 



where x{t) is a solution of (|i|) for a given realization of s and for a given initial condition. The average is taken over 
the noise realizations. Taking the time derivative of (A_2) and using eq. ([T]) one gets 

(A3) 



dtQix,t) = -9, [v,{x,t)Qix,t)] + d,{s,{t)S{x{t) - x)). 



Taking advantage of the Gaussian nature of s the average in (A3) can be rewritten as 



^t&{x,t)^-^^[v^{x,t)Q{x,t)]+ S 



dt' 



to 



S{x{t) — x] 



(A4) 



Because of the (5-function a closed equation is possible only if the functional derivative either does not involve the 
process x or depends on it solely at the "Markovian" end-point t = t' . At this stage thus we cannot simplify the 
master equation any further. In the limit of small correlation time r a closed equation can be derived by performing 
the largely used small r expansion. If the noise is close to the white noise limit (r = 0) it is reasonable to expand 
the functional derivative about its Markovian value, i.e., the one obtained for the (5-correlated noise. The Taylor 
expansion of 5xj{t) / 5si{t') around the Markovian end point t' = t is 



dxj{t) ^Xj{t) 



d Sxj (t) 



dt' Ss,{t') 



(t'-t) + 



^S,j-d,vj{x,t) (t'-t) + 



(A5) 
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Inserting the expansion (A5) into (A4), keeping only the first terms in r and negletting the transients, i.e. the letting 
to — oo we obtain after a straightforward algebra the small t master equation: 



where 



dtQ{x,t) = -5, [v,{x,t) e{x,t)] + df. [V,j{x,t) Q{x,t)] 



'Dij{x,t) = % + ^ {diVj{x,t) + djVi{x,t)) 



(A6) 



(A7) 



and use of incompressibility has been made. 

We note that this expansion does not converge uniformly in x, and the diffusion coefficien t (lATt ) may exhibit 
negative values, thereby introducing unphysical singularities into the problem. In other words ( |A6| ) and (A7) does 
not constitute a truly Markovian process with well-defined corresponding Langevin equation driven by white noise. 
In general these equations are valid only for r ^ 1 and t/ct^ <C 1 (in dimensionless units). 



APPENDIX B: THE GENERALIZED FORMULAS IN THE STEADY SHEAR FLOW CASE 



The exact expression of the eddy diffusivity tensor for a stationary bidimensional shear flow 

V = {v{y), 0) 

with v{y) a periodic function in y can be deduced from ( |3^ ) and reads 

1 



D22 = cr . 



(Bl) 



(B2) 



Applying the generalized formulas (| 
For the shear flow eq. (BSl) reads 



) and (B^) to the shear flow we want to obtain the expressions (B2) 



{dt + vdi~a'^d'^)w^°^ = -v 
{dt + vdi- a^d^)w^^^ = 0. 



(B3) 
(B4) 



Inserting the stationary solution of eq. (B4), w'^'' = 0, into eq. (B4) and taking the Fourier transform v{k) of v{y) 
we obtain: 



0-2 /c2 



Thus eq. (|66|) becomes 



{dt +vdi~ = (p _ q)a^dlv 

{dt + vdi- (t292)u;W ^ 0. 



The stationary solutions are 



w[^\y) = -{p-i) v{y) 



= 0. 



Using eq. the eddy diffusivity tensor is 
Dg = a' 



and 



T{p~2q)l-Jdk \v{k) p+0(r2). 



(B5) 

(B6) 
(B7) 

(B8) 



(B9) 
(BIO) 

(Bll) 



(B12) 
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A comparison between (B12) and (B2) shows that the generaUzed equations lead to the exact expression of the eddy 



diffusivity tensor for a stationary shear flow if 

p = 2q+l. 

The same condition is found if we consider a time-dependent shear flow. 
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FIG. 1. The ratio A as a function of the cr^ for the 
two-dimensional AB flow with A = B = 1. The continuous 
line is the prediction obtained by while the points are the 
numerical results obtained from the (^). All the quantities 
here are supposed to be dimensionless. 
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FIG. 2. The streamlines for the quasi two-dimensional 
flow with k = 2. 
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FIG. 3. The streamlines for the two-dimensional AB flow 
( [Fil ) with A = B — 1. X, y, A and B are supposed to be all 
dimensionless. 
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